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LIQUID FUEL SPRAY PROCESSES IN HIGH-PRESSURE GAS FLOW 


Robert D. Ingebo* 


Atomization of single liquid jets injected down- 
stream in high-pressure and high-velocity airflow 
was investigated to determine the effect of air- 
stream pressure on mean drop size as measured 
with a scanning radiometer developed at NASA 
Lewis Research Center. For aerodynamic-wave 
breakup of liquid jets # it was found that mean 
drop diameter D 32 varied with airstream pres- 
sure P a as follows: D 32 ' • This 

result was compared with the effect of airstream 
pressure on drag and heat-transfer coefficients 
of accelerating and vaporizing diops. 


INTRODUCTION 


The fluid mechanics of producing sprays and the aerothermodynamics 
of droplet acceleration and vaporization have been receiving con- 
siderable attention in spray process research studies. Spray 
drying for food processing, spray painting of manufactured items, 
pesticide sprays for agricultural use, icing sprays in wind tun- 
nels, and fuel sprays in gas turbine combustors are some of the 
major fields of application of spray process studies. In the 
field of combustion, Ingebo (1) determined that the fuel spray 
vaporization rate is very important in determining burning 
characteristics such as ignition delay, flame stability, and 
completeness of burning. Fuel droplets formed by atomization 
simultaneously accelerate and evaporate in airstreams and thereby 
produce a combustible ruel-air mixture. Chemical reaction between 
fuel and air, or oxidant vapor, produces a combustion gas stream 
of relatively high velocity, temperature, and pressure, which 
accelerates the drops as the chemical reaction continues to accel- 
erate the combustion gas flow in the combustor. Thus, a clear 
understanding of liquid fuel spray combustion requires knowledge 
of how aerothermodynamics and fluid mechanics may control the 
processes of atomization, vaporization, and acceleration, and how 
these processes are related to chemical reaction in determining 
combustor performance and exhaust emissions. More specifically, 
initial drop size and velc-city as well as combustion gas condi- 
tions must be known to calculate heat-transfer and drag coeffi- 
cients and thereby determine the lifetime histories of burning 
fuel drops. New data of this type are especially needed to 
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determine rates of drop vaporization and combustion under high- 
pressure gar flow conditions. 

Empirical expressions are quite useful in the design of fuel 
injectors for combustor applications since fuel nozzles can be 
calibrated with water, and mean drop diameters can then be calcu- 
lated for the desired fuel, such as Jet-A, with the aid of such 
correlations. The product of the Weber and Reynolds numbers 
WeRe has been correlated with atomization data, and, for the 
condition WeRe > 10°, atomization occurs in the regime of aero- 
dynamic breakup. This condition was chosen for this investigation 
since the airstream velocity in turbojet combustors at idle, take- 
off, and cruise operating conditions is usually quite high; i.e., 
combustor reference velocities are generally of the order of 20 to 
30 m/s. 

Since previous studies on atomization were limited to atmos- 
pheric pressure conditions, the present investigation was made to 
determine the effect of airstream pressure on mean drop size. 

Water jets injected downstream were atomized in the ai r< "ow, and 
mean drop diameters were measured with an improved scanning radio- 
meter recently developed at NASA Lewis Research Center. Combustor 
inlet-air static pressure was varied from 0.10 to 2.10 tJPa over a 
range of airflow rates per unit area of 10 to 177 g/(cm-} (s) at 
293 K. Two injector tubes with orifice diameters of 0.1016 and 
0.216 cm, respectively, were used with liquid flow rates of 27 to 
68 liters/hr. Mean drop diameter data were correlated with the 
Weber-Reynolds number WeRe and a pressure sensitive molecular 
scale dimensionless group gl/c 2 . The resulting empirical rela- 
tions for liquid jet breakup were compared with expressions 
obtained in previous experimental investigations. Also, the 
effect of airstream pressure on atomization is compared with air- 
stream pressure effects on droplet acceleration and vaporization 
rates as reported in previous spray studies in order to bring 
together high-pressure airflow data for spray processes reported 
in the literature. 


APPARATUS AND PROCEDURE 

The high-pressure airflow apparatus and auxiliary equipment are 
shown in figure 1. The test section, water spray, and scanning 
radiometer optical path are shown in figure 2. Two 5.1-cm- 
diameter windows pass the laser beam through the spray 25.4 cm 
downstream of the fuel injector. Single jets of water were 
injected downstream at the duct centerline with the fuel tube 
shown in figure 2. Two different 0 . 635-cm-outside-diameter tubes 
were used to obtain drop size data with two orifice diameters of 
0.1016 and 0.216 cm, respectively. 

At 273 K, single water jets were injected downstream at flow 
rates of 27 and 68 liters/hr, respectively, as measured with a 
turbine flowmeter. After setting water and airflow rates, mean 
drop diameter data were obtained with the scanning radiometer 
shown in figure 2. The optical system consisted of a 1-mW helium- 
neon laser, a 0.003-cm-diameter aperture, a 7 . 5-cm-diameter colli- 
mating lens, a 10-cm-diameter converging lens, a 5-cm-diameter 
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collecting lens, a scanning disk with a 0 . 05-cm-diameter hole, a 
timing light, and a photomultiplier detector. 

The spatial resolution of the scanning radiometer is 2.86 cm, 
which is also the laser beam diameter in the test section where a 
sufficient volume of spray was sampled to minimize spray pattern 
effects on mean drop diameter. Also, the effect of drop size 
distribution function on the scanning radiometer measurements is 
discussed in detail by Dobbins, et al. (2) , who found that the 
irradiance distribution is only weakly related to the particle 
diameter distribution function. Thus, irradiance distribution 
was used in determining Sauter mean diameter, and changes in the 
drop size distribution function were assumed to have negligible 
effect on mean drop size correlations. Reproducibility tests 
showed that experimental measurements of mean drop diameter agreed 
within +5 percent. Three sets of monosized polystyrene spheres 
having diameters of 25, 50, and 100 pm, respectively, were used 
to calibrate the scanning radiometer. A more complete description 
of it and the method of determining mean particle size are 
discussed by Buchele (3,4). 


EXPERIMENTAL RESULTS 

To determine the effect of airstream pressure on liquid jet atomi- 
zation, an improved model of the scanning radiometer was developed 
at NASA Lewis Research Center. It was used to measure the mean 
diameter of droplets in sprays produced by liquid jets breaking up 
in high-pressure and high-velocity airflow. 

Liquid jet atomization 

In a previous study conducted with atmospheric pressure air- 
streams, drop size distribution data for the breakup of liquid 
jets were obtained with the high-speed droplet tracking camera 
shown in figure 3 and described by Ingebo and Foster (5) . Photo- 
micrographs of the droplets and a schematic of the liquid jet 
breakup process are shown in figure 4. Short wavelength distur- 
bances near the injector orifice produce relatively small droplets 
and grow into longer wavelengths as the jet penetrates farther 
into the airstream and produces relatively large drops. Mean 
drop size measurements obtained with this imaging technique gave 
the data shown in figure 5 for five different liquids. The fol- 
lowing correlation of the ratio of mean diameter to orifice diam- 
eter D3o/ D o with the Weber-Reynolds number WeRe was derived 
for values of WeRe < 10°: 

d 30/ d o = 3.9 (WeRe) _0 * 25 (1) 

In deriving this relationship in reference (5) , the spatial size 
distribution data obtained with the high-speed camera were mass 
weighted using the liquid-concentration profile data obtained 
with a sampling probe. The temporal size distribution produced 
is more useful in combustion studies. 

In a later study, an early model of the scanning radiometer as 
described by Ingebo (6) was used to study liquid jet breakup. 

Here it was found that equation (1) was valid only in the regime 
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of capillary-wave breakup which occurs for values of WeRe < 10*>. 
At values of WeRe > 10 , it was found in reference (6) that 

Do/D 32 = 0.027 (WeRe) 0 - 4 (2) 

for liquid jet atomization in the regime of acceleration-wave 
breakup. Date from reference (6) for both capillary and 
acceleration-wave breakup are plotted in figure 6. 


In this investigation an attempt was made to extend the appli- 
cation of equation (2) to the atomization of liquid jets in high- 
pressure airflow. Single jets of water were injected axially 
downstream in airstreams having static pressures from 0.1 to 
2.7 MPa. Atomization data obtained with the two different fuel 
injector tubes are shown in figure 7 where the reciprocal mean 


drop diameter D, 


-1 


m 


is plotted against airstream mass velocity 


PaV r 


on 




1 decreases mark- 


PaV r . As shown, the effect of 

ealy as airstream pressure is increased from 0.1 to 2.7 MPa. To 
extend equation (2) and make it applicable to high airstream pres- 
su'es, (Dq/Djp) (WeRe) * 4 is plotted against the pressure sensitive 
grt up gl/c^(fig. 8). In this dimensionless group g is the 
acceleration due to gravity, 1 the mean free molecular path, 


and 


the root-mean-square velocity of the air molecules. 


The evaluation of this dimensionless group is as follows. The 
mean free molecular path 1 may be expressed as 

1 = l/v'InDg 2 = 6.11xl0 -6 cm 

since the number of molecules per cubic centimeter n is 2.7x10^-^ 
at a temperature of 0 °C and a pressure cf 1 atm and the diameter 
of an air molecule Dg is 3.7x10“^ cm. 

The root-mean-square molecular velocity c may be expressed as 
c = ( 3RT/M) 0 • ^ , w hich yields 


5 2 . 13H_ 8.31xl0 7 . ) . ( 2 73) . 2. 35xl0 9 (cm/s) 2 

Since g is 980 cm/s/s, 

gl/S 2 = (980) (6.11x10-**) = J-55xl0 -12 
2.35xlO y 

From the data plotted in figure 8, it is evident that 
D m “ 1 % (gl/c^) Thus, the following expression is obtained 
from the data plotted in figure 9: 


D o /D 32 = l-2(WeRe)°- 4 (gl/5 2 ) 0 * 15 


( 3 ) 


which agrees with the expression obtained by Ingebo (7) for the 
cross stream injection of liquid jets - i.e., 

D 0 /D 32 « 1.4 (WeRe) 0 * 4 (gl/c 2 ) 0 ‘ 15 
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for acceleration-wave breakup of liquid jets with values of 
WeRe > 10 . Since the molecular-scale momentum-transfer group 
gl/c 2 is inversely proportional to airstream pressure, it is 
evident that gl/c 2 -v P a "^ . It was just shown that, at 273 K 
and 1 atm pressure, gl/c 2 = 2.55xl0“* 2 . Also, since WeRe ^ P a , 
equation (3) shows that D 32 'v. P a ” 0,2 . Thus, the effect of 
airstream pressure on D 32 is somewhat less than that which 
would be obtained by extrapolating an atmospheric pressure expres- 
sion up to high airstream pressures. For example, equation (2) 
would indicate that D 32 P a “° . 

A comparison of the results obtained by several investigators 
of airstream effects on atomization is given in table I. Best 
agreement was obtained with Lorenzetto and Lefebvre ( 8 ) where it 
was found that D 32 'v P a '®‘ . This is in good agreement with 
equation (3) since it is assumed that D m = D 32 for the measure- 
ments of D_ made with the improved scanning radiometer in the 
present stuay. Kim and Marshall (11) and Wolfe and Andersen (12) 
give values more in agreement with the theoretically predicted 
relationship D 32 x P a ” 0,67 given by Adelberg (13). The atomi- 
zation data presented in this study and in reference ( 8 ) were 
obtained over a much broader airstream pressure range than that 
used in references (11) and (12). This indicates a need to 
improve the theory proposed in reference (13) , which does not 
account for molecular-scale forces such as those given in the 
group gl/c 2 . 

Droplet drag coefficients for accelerating sprays 

The high-speed droplet tracking camera used to obtain drop size 
data also gave drop velocity data. Thus, it was possible to 
determine drag coefficients of drops accelerating in gas streams. 
In the study described by Ingebo (14) , it was found that the drag 
coefficient C,j could be correlated with Reynolds Re ased on 
drop diameter as follows: 

Cd = 27Re -0 * 84 


at atmospheric pressure. 

Studies of propellant atomization and droplet acceleration in 
the hiqh-pressure rocket combustor shown in figure 10 were made 
with the droplet tracking camera as described by Ingebo (15) . It 
was found that the drag coefficient based on the mean drop diam- 
eter D 20 could be correlated with the Reynolds number Re and 
the dimensionless group gl/c^ as follows: 


Cd 


0 . 13Re -0 * 84 (gl/c 2 ) ~ 0 ’ 2 


(4) 


where Re is based on the area number mean drop diameter £> 20 * 

A plot of the data used to derive this expression is shown in 
figure 11. From equation 4, the effect of airstream pressure on 
drag coefficient appears to be C<f ^ P a ’ 0,6 • As in the study 
of atomization, the effect of airstream pressure is somewhat less 
than that predicted from atmospheric pressure studies. 
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Heat-transfer coefficients for vaporizing drops 


In the study of vaporization rates of single drops in atmos- 
pheric pressure airstreams, Ingebo (16) found that the Nusselt 
number Nu could be correlated with the product of the Reynolds 
and Schmidt numbers ReSc as follows: Nu = 2 + 0.30 (ReSc) 

where Nu = hD/kg, h is the heat-transfer coefficient, D the 
drop diameter. Sc = u/b gw , k g and p the gas thermal conductiv- 
ity and viscosity, respectively, and b gv the vapor diffusivity. 

In a later study by Ingebo (17) , at an airstream static pressure 
range of 0.6 to 2 atm, data for the Nusselt number were obtained as 
shown in figure 12 and the following expression was derived: 

Nu = 2 + 2.6xl0 6 [(Re&c) (gl/c 2 )] 0 , 6 (k g /k v ) 0,5 (5) 

where k v is the thermal conductivity of the vapor. Thus, 

Nu 'v p a 0. However, airstream pressure still influenced vapori- 
zation rate by its effect on temperature difference across the 
air-vapor film surrounding the drop. 


CONCLUSIONS 


As a result of this study and previous spray studies in high- 
pressure airstreams, graphical or numerical integration methods 
of calculating spray vaporization rates will not be very useful 
until more atomization, drag coefficient, and heat-transfer 
coefficient data are available for very high-pressure airstreams. 
The effects of airstream pressure on mean drop diameter D 32 , 
drag coefficient C^, and Nusselt number Nu are shown in 
table II for comparison. Probably the most critical need at pres- 
ent is for high-pressure airstream data on drop size distribution 
and on mean drop size for spray processes in high-temperature and 
highly turbulent gas streams. Also, a more fundamental picture 
of simultaneous heat, momentum, and mass transport properties is 
needed for a more comprehensive treatment of fuel spray processes 
and the overall combustion process. 


NOMENCLATURE 

bg, w = molecular mass diffusivity (g/(cm)(s)) 

Cd = drag coefficient 

c = root-mean-square molecular velocity (cm/s) 

D = diameter (cm) 

D 32 = Sauter mean diameter, ^nD-VjnD^ (cm) 

g = acceleration due to gravity (980 cm/s^) 
h = heat-transfer coefficient (g-cal) / (sec) (cm 2 ) (°C) 
k = thermal conductivity (g-cal) / (sec) (cm^) (°C/cm) 

1 = mean free molecular path (cm) 
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M = molecular weight 

Nu = Nusselt number for heat transfer, hD/k 
n = number of molecules per unit volume 
P = static pressure (atm) 

R = universal gas constant (8.31rl0 7 ergs/ (K) (mole) ) 

Re = Reynolds number based on orifice diameter, D 0 v c /y^ 

Sc = Schmidt number based on mass diffusivity, v/bg,w 
V = velocity (cm/s) 

We = Weber number based on orifice diameter, D Q P a V r /a 
T = kinematic viscosity (cm 2 /s) 

a = surface tension (dynes/cm) 

p = density (g/cm^) 

y = absolute viscosity (g/(cm)(s)) 

Subscripts : 
a = airstream 

g = gas molecule 

1 = liquid 

m = mean 

o * orifice 

r = relative 

v = vapor 
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TABLE I. - EXPERIMENTALLY DETERMINED E) ’ONENTS OF 

PROPERTIES PROPORTIONAL TO RECIPROCAL MEAN DROP 

DIAMETER D* 1 FOR AERODYNAMIC -WAVE BREAKUP 
m 

OF LIQUID SHEETS AND JETS 


Sources for 
exponents 

Airstream 
velocity , 

v a/ 

cm/s 

Orifice 
diameter , 
Do • 
cm 

Airstream 

pressure, 

9 

MPa 

Experimental (eq. (3)) 

1.2 

0.20 

0.25 

Nukiyama and 
Tanasawa (9) 

1.0 

.00 

0 

Weiss and Worsham (10) 

1.33 

-.16 

(a) 

Kim and Marshall (11) 

1.14 

.00 

.57 

Wolfe and Andersen (12) 

1.33 

b - . 17 

.67 

Lorenzetto and 
Lefebvre (8) 

1.0 

0 i 

1 

f 

.30 


a Used variable 1 + p a /pi. 

^Exponent for initial drop diameter instead of orifice 
diameter. 


TABLE II. - COMPARISON OF EFFECT OF AIRSTREAM 
PRESSURE P„ ON SPRAY PROCESS PARAMETERS 
D m and C d AND NUSSELT NUMBER Nu 


Atomization 

D /D = 1.2(WeRe)°- 4 (gl/5 2 ) 0 ' 15 
-0 25 

D x P 
m a 

Acceleration 

1 

c d = 0.:3Re‘°* 84 (gl/c 2 )" 0,2 

C, -v p -0 ’6 4 
d a 

Vaporization 

Nu = 2 + 2.6xl0 6 [ (ReSc) (gl/c 2 )] 0,6 

h -v, P° 
a 
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Figure 5. - Relation between mean to orifice diameter ratio 
and Weber-Reynolds number ratio. 
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Figure 6. - Relation between orifice to mean diameter ratio and Weber-Reynolds 
number. 
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Figure 7. - Variation of reciprocal mean I 

drop diameter with airstream mass 
velocity. 
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(b) Orifice diameter, D Q , 0. 1016 cm. 

Figure 8. - Correlation of mean drop diameter D m 
with airstream pressure sensitive group gl/c - 2L 
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Figure 12. - Correlation of Nusselt number with momentum-transfer groups 
and thermal-conductivity ratio. 



